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ABSTRACT 

The existence of X-ray luminous gaseous coronae around massive disc galaxies is a long- 
standing prediction of galaxy formation theory in the cold dark matter cosmogony. This pre- 
diction has garnered little observational support, with non-detections commonplace and de- 
tections for only a relatively small number of galaxies which are much less luminous than 
expected. We investigate the coronal properties of a large sample of bright, disc-dominated 
galaxies extracted from the GIMIC suite of cosmological hydrodynamic simulations recently 
presented by Crain et al. Remarkably, the simulations reproduce the observed scalings of X- 
ray luminosity with if -band luminosity and star formation rate and, when account is taken of 
the density structure of the halo, with disc rotation velocity as well. Most of the star formation 
in the simulated galaxies (which have realistic stellar mass fractions) is fuelled by gas cooling 
from a quasi-hydrostatic hot corona. However, these coronae are more diffuse, and of a lower 
luminosity, than predicted by the analytic models of White & Frenk because of a substantial 
increase in entropy at z ~ 1 — 3. Both the removal of low entropy gas by star formation and 
energy injection from supernovae contribute to this increase in entropy, but the latter is dom- 
inant for halo masses Af 2 oo ^ 10 125 M Q . Only a small fraction of the mass of the hot gas 
is outflowing as a wind but, because of its high density and metallicity, it contributes dispro- 
portionally to the X-ray emission. The bulk of the X-ray emission, however, comes from the 
diffuse quasi-hydrostatic corona which supplies the fuel for ongoing star formation in discs 
today. Future deep X-ray observations with high spectral resolution (e.g. with NeXT/ASTRO- 
H or IXO) should be able to map the velocity structure of the hot gas and test this fundamental 
prediction of current galaxy formation theory. 
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1 INTRODUCTION 

The rarity of X-ray detections of hot gaseous coronae surrounding 
disc galaxies poses a fundamental challenge to the current view of 
how galaxies form. Whilst the possibility that galaxies might con- 
tain extra-planar gas was already raised by Spitzer (1956), it was 
White & Rees ( 1978) who first proposed that hot gaseous coronae 
are an integral part of the galaxy formation process. Their two-stage 
theory posited that hot gas reservoirs build up as gas condensing 
onto dark matter haloes is heated by thermodynamic shocks and 
adiabatic compression. Subsequent radiative cooling of this gas es- 
tablishes a 'cooling flow' that fuels ongoing star formation. 

The ubiquity of extended soft X-ray emission from the cool- 
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ing of galactic gaseous coronae was first predicted by |White &] 
|Frenk| ( fT991| hereafter WF91), whose model built upon the ideas 
sketched by White & Rees and established the foundation for inter- 
preting galaxy evolution within the cold dark matter cosmogony. In 
their framework, the dark matter haloes that host present-day disc- 
dominated L* galaxies engender gas density profiles that are con- 
ducive to efficient radiative cooling, via line-emission and thermal 
Bremsstrahlung, out to radii beyond the optical extent of the central 
galaxy. The associated cooling rate is sufficient to fuel ongoing star 
formation in disc galaxies, and thus offset the disruption of discs 
by mergers (Walker et al. 1996, Barnes 
( pfstathiou et al.|1982||Mo et al.|1998 



1998) 



and bar-instabilities 



Syer et al. 1999), enabling 



the model to reproduce the abundance of morphological types ob- 
served in the local Universe (e.g. Driver et al. 2006 . 2007 , Fukugita 



|et al.|2007||Bernardi et al.|2010[|Parry et al.|2009 
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The radiation associated with cooling from gas in the haloes 
of L* galaxies is predicted to fall primarily in the soft X-ray band 
and to have typical surface brightnessess that are readily observ- 
able with the XMM-Newton and Chandra telescopes. This predic- 
tion is compelling because analytic and semi-analytic models based 
upon the WF91 framework successfully reproduce a broad range 
of galaxy properties, such as the luminosity function in the opti- 
cal ( |Cole et al,Tl9 94; Kauffmann et al. 1999a, Some rville & Pri-| 
mack 1999; Cole et al. 2000}, infrared jLace y et al. 2008) and sub- 
millimetre (Baugh et al. 20051 wavebands; the bimodality in the 
colour-magnitude plane; the apparent 'downsizing' of galaxy for- 
mation; and the black hole scaling relations jBower et al.||2006| 
|Croton et al.|2006a||De~Lucia et al.|2006) . 

The existence of X-ray luminous coronae around massive el- 



liptical galaxies is well established (e.g. Forman et al.|1985| Kim 



eTaTl[l992l |Q' Sullivan eraLl[200T] |David et al.||2006l |Sun et~ 
2007||Jeltema et a l. 2008, Sun et al.|20"09") . However, this is gener- 



ally not regarded as a test of the WF91 model because this difffuse 
X-ray emission could be associated with gas returned to the ISM 
through stellar evolution (i.e. supernovae and AGB stars, e.g. [Read] 
|& Ponman|1998||Mafhews & Brighenti|1998") or with hot gas con- 
fined by the potential of galaxy groups - the most common envi- 
ronments of massive ellipticals (e.g. Ponman et al. 1994, Trinchieri 
|et al.|1994|[Mulchaey et al.|1996||Trrn"chieri et al.|1997) . The search 
for diffuse soft X-ray emission from hot coronae around massive, 
isolated disc galaxies is therefore the most direct test of the canon- 
ical galaxy formation picture. 

Adopting this premise, |Benson et al.| ( |2000| hereafter BOO) 
analysed ROSAT observations of three luminous, nearby galaxies, 
of which two (NGC 2841, NGC 5529) are disc-dominated, whilst 
the third (NGC 4594, the 'Sombrero galaxy') is a bulge-dominated 
Sa. They found no convincing evidence for diffuse X-ray emission: 
their upper limits on the soft X-ray luminosity are over an order 
of magnitude below the luminosity predicted by the WF91 cooling 
flow model. Similarly, an updated analysis of Chandra observa- 
tions of NGC 5746 by Rasmussen et al. (2009 1 failed to produce a 
significant detection. The greater sensitivity of XMM-Newton and 
Chandra relative to ROSAT has yielded the detection of a diffuse 
X-ray component for NGC 4594 JLi et al.|2006l|2007] l, but even in 
this case the luminosity shortfall relative to the cooling flow model 
is approximately two orders of magnitude. Moreover, energy in- 
put from supernovae has been claimed as the most likely source 
of the emission in this case. Detections of diffuse emission from 
disc-dominated galaxies are, however, now becoming more com- 
mon (e.g.|Strickland et al.|2004||Wang|2005||TTin"mann et al.|2006] 



|Owen & Warwick|2009| |Sun et al.|2009| > and, in general, exhibit 
a correlation with the star formation rate. A small number of spec- 
tacular cases exhibit X-ray emission that has an obvious biconi- 
cal morphology driven by nuclear starbursts. This has led to the 
common interpretation that heating by supernovae is the dominant 
(and perhaps sole) mechanism by which local disc galaxies gener- 
ate extra-planar X-ray emission, rather than the thermalisation of 
gravitational potential energy. 

Another potential complication is the emission originating 
from point sources (X-ray binaries in particular), which can con- 
tribute a significant fraction of the total X-ray flux. However, this 
is much less of a problem for the current generation of X-ray satel- 
lites, which have excellent angular resolution (Chandra in particu- 
lar) that enables identification and masking of bright point sources, 
and high spectral resolution, which allows the removal of the con- 
tribution of unresolved, faint point sources (but see |Revnivtsev| 
|et al.|2008). Typically, the removal of unresolved point sources is 



achieved by including a power-law component, in addition to the 
thermal component, when fitting the X-ray spectrum. 

The dearth of convincing X-ray detections of hot coronae is 
puzzling and suggests that the canonical view of disc galaxy for- 
mation is incomplete. Analytic models such as WF91 necessar- 
ily rely on simplification of the complex and non-linear processes 
that are the essential elements of galaxy formation, such as metal- 
dependent gas cooling, anisotropic gas dynamics, star formation, 
feedback and chemical enrichment. Potentially, a more general 
view of the formation of hot coronae and the growth of galaxies 
is offered by direct, hydrodynamical simulations. 

|Toft et aT| ( |2002[ > highlighted a discord between hydrody- 
namic simulations and analytic models, finding that efficient ra- 
diative cooling can lead to the mass of hot gas confined by dark 
matter haloes being substantially lower than the predictions of cool- 
ing flow models. Since X-ray emissivity is sensitive to gas density 
(oc p 2 ), the overall X-ray luminosity of the galaxy is commensu- 
rately reduced; their model produced galactic coronae with soft X- 
ray luminosities two orders of magnitude lower than predicted by 
simple cooling flow models, and consistent with the limited obser- 
vations available at the time. Interestingly, this result is in contrast 
to the naive expectation of analytic models where, in general, it is 
assumed that cooling efficiency and X-ray luminosity are closely 
correlated. 

It is worth noting that Rasmuss en et al.| (J2009 ) recomputed 
the evolution of the systems presented by Toft et al. (20021, us- 
ing an updated simulation code and higher-resolution initial condi- 
tions. The more recent study reported X-ray luminosities that were 
a factor of ~ 2 lower than found by Toft et al. (2002). Whilst these 
studies highlight the value of hydrodynamic simulations as a means 
to test the validity of the simplifications made in analytic models, 
their differences point to the fact that predictions derived from sim- 
ulations are also subject to significant uncertainties stemming from 
resolution effects and the algorithmic implementation of the under- 
lying physics. 

In this paper, we investigate the X-ray halo problem using the 
very large Galaxies-Intergalactic Medium Interaction Calculation 
(GIMIC; |Crain et al.|2009| hereafter C09), a suite of high resolution 
hydrodynamic resimulations of regions drawn from the Millennium 
simulation (Springel et al. 2005). A central aim of this project is to 
study the interaction of galaxies with their gas haloes and with the 
external intergalactic medium (IGM). The GIMIC simulations are 
well suited to our purposes, not only because they feature many 
well-resolved L* galaxies with stellar discs, but also because their 
convergence behaviour is well understood (see C09), and because 
the only parameter that was tuned is the mass loading of winds, 
chosen to ensure a reasonable match to the observed cosmic star 
formation rate density. 

The highest resolution realisations of the simulations repre- 
sent the stellar discs of galaxies similar to the Milky Way with 
~ 10 5 particles, whilst the intermediate resolution realisations do 
so with a factor of 8 fewer. In both cases, the resolution is suf- 
ficient to classify the morphology of the galaxies as being either 
disc- or spheroid-dominated, enabling us to subsample the simu- 
lated galaxy population so as to most closely match observational 
samples. Combined, the intermediate-resolution simulations fol- 
low a comoving volume of 1.6 x 10 5 (h -1 Mpc) 3 , yielding ap- 
proximately 460 galaxies at z = with stellar mass greater than 
10 10 /i -1 M©. The simulations include mass and energy feedback 
due to supernovae, and chemical enrichment from type la and type 
II supernovae and asymptotic giant branch (AGB) stars. Radiative 
cooling is computed on an element-by-element basis under equilib- 
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rium conditions in the presence of a realistic photoionising UV/X- 
ray background. This allows us to investigate whether X-ray emis- 
sion is generated primarily by the presence of a gravitationally sup- 
ported coronae, or by extra-planar gas heated by supernovae. We 
are also able to assess the importance of metal cooling to the emis- 
sion. 

This work is part of the programme of the Virgo consortium 
for cosmological simulations. The simulations adopt the same cos- 
mological parameters as the Millennium simulation (Springel et al. 
[20Q5] >: f! m = 0.25, fi A = 0.75, fi b = 0.045, n s = 1, as = 0.9, 
Ho = 100 /ikms" 1 Mpc^ 1 , h — 0.73. Our choice of cosmolog- 
ical parameters is consistent with the most recent determinations 
from the cosmic microwave background radiation and other cos- 
mological tests, with the exception of erg, which is higher by ~ 2a. 
We do not expect this difference to affect the results and conclu- 
sions of this work, since for the halo mass range considered here it 
results in only a small change in the epoch of formation. 

The paper is laid out as follows. In §[2] we describe our meth- 
ods, giving a brief overview of the GIMIC simulations and our main 
postprocessing techniques. In §[3] we confront the simulations with 
observational constraints, whilst in §p]we explore the analytic the- 
ory and compare it with the results of the hydrodynamic simula- 
tions. We summarise and discuss the results in §|5] A short appendix 
is included, in which we discuss the sensitivity of our numerical re- 
sults to resolution. 



2 METHODS 

2.1 Simulations 

The GIMIC simulations are described in detail in C09, where thor- 
ough discussions of the generation of the initial conditions, the sim- 
ulation code, and the initial results may be found. Here, we present 
only a brief overview and limit the description to aspects that are 
specifically relevant to this study. 

GIMIC is designed to circumvent the unfeasibility - due to 
computational expense - of simulating large cosmological volumes 
(L > 100 ft -1 Mpc) at high resolution (m gas < 1O 7 /i _1 M ) 
to z — 0, Using 'zoomed' initial conditions ( |Frenk et al.||1996| 
|Power et al.|2003| |Navarro et al.|2004| >, GIMIC follows, with full 
gas dynamics, the evolution of five roughly spherical regions drawn 
from the Millennium simulation. In order to encompass a wide 
range of large-scale environments, the regions were chosen such 
that their overdensities deviate by (—2, — 1, 0, +l,+2)cr from the 
cosmic mean, where a is the rms mass fluctuation, on a scale of 
18 h" 1 Mpc, at z — 1.5. The +2a region was additionally con- 
strained by the requirement that it be centred on a rich galaxy clus- 
ter halo. In practice, this ensures that the simulations include rare 
cosmological features, since the — 2a region is also approximately 
centred on a sparse void. Each region has an approximate comoving 
radius of 18 h^ 1 Mpc, except the +2a region which was enlarged 
to a radius of 25 h^ 1 Mpc in order to accommodate the rich clus- 
ter. The remainder of the 500 3 (ft -1 Mpc) 3 Millennium simulation 
volume is modelled with collisionless particles at much lower res- 
olution in order to follow the large-scale structure. The scheme is 
illustrated schematically by Fig. 1 of C09. 

The initial conditions were realised at intermediate resolu- 
tion (m gas = 1.16 x 10 7 /i" 1 Mq) and high resolution (m gas = 
1.45 x 10 6 ft -1 Mq); we reserve the term 'low resolution' for the 
original Millennium simulation, in which the collissionless parti- 
cles, representing a composite of baryonic and dark matter, have 



mass 8.6 x 10 8 h~ x Mq. Gravitational forces on the baryonic and 
high-resolution dark matter particles were softened over an identi- 
cal length scale; we adopt a softening length that is initially fixed 
in comoving coordinates, but becomes fixed in physical coordi- 
nates at a predefined redshift, i.e. e CO m(a) = min(e com , e™j^f s /a). 
The softenings were chosen such that at z — 3, they are fixed 
at e ™h^s = (1.0,0.5) ft _1 kpc for the intermediate- and high- 
resolution runs, respectively. 

The simulations were performed with the TreePM-SPH code 
GADGET3, a substantial upgrade of GADGET2 (Springel 2005} that 
includes: 

(i) a new domain decomposition algorithm that improves load 
balancing, particularly for simulations with strongly clustered par- 
ticle distributions run on parallel supercomputers with a large num- 
ber of cores (Springel et al. 2008); 

(ii) a prescription for star formation designed to enforce a local 
Kennicutt-Schmidt law ( |Schaye & Dalla Vecchia|2008) ; 

(iii) the contribution of metals to the cooling of gas, computed 
element-by-element, in the presence of an imposed UV-background 
( |Wiersma et al.|2009a| >; 

(iv) stellar evolution and the associated delayed release of 11 
chemical elements (Wiersma et al. 2009b); 

(v) galactic winds that pollute the IGM with metals and can 
quench star formation in low-mass haloes (Dalla Vecchia & Schaye 
12008) 1. 

The code does not, however, model the evolution of black 
holes or feedback effects associated with them. The hydrodynamics 
implementation, taken from GADGET2, is the entropy conserving 
formulation of smoothed particle hydrodynamics (SPH; Gingold & 
Monaghan 1977||Lucy|1977) , as discussed in Springel & Hernquist 
(2002). For further details, the reader is referred to C09 and refer- 
ences therein. All five regions were evolved to z = at intermedi- 
ate resolution. Only the —2a region was evolved to z = at high 
resolution (the — la, 0a and +ler regions were run to z — 2; the 
+2cr region was not run at high resolution). Since we are interested 
in the low redshift evolution of galaxies, we follow the philosophy 
of C09 and consider the five intermediate-resolution runs as a fidu- 
cial simulation set, and use the high-resolution realisation of the 
—2a simulation to assess how sensitive our results are to numeri- 
cal resolution. An assessment of numerical convergence is supplied 
in the Appendix. 

Owing to its unique initial conditions, GIMIC represents a 
complementary simulation suite to the Overwhelmingly Large 
Simulations (OWLS; |Schaye et al.|201 0). Both projects are based 
on the same simulation code. The aim of OWLS is to investigate the 
dependence of various properties on the parametrisation of subgrid 
physics, whilst GIMIC accesses a wider range of environments for 
a single physics implementation of OWLS (the MILL simulation in 
OWLS). 

We note that we are concerned here with the establishment 
of hot coronae, for which shock heating (induced by gravitational 
collapse and the impact of galactic winds with the intergalactic 
medium) can play a key role. The ability of our hydrodynamics 
scheme to model this process robustly is therefore important. SPH 
has been shown to reproduce the analytic solutions of simple prob- 
lems such as shock tubes, spherical collapse, and Sedov blasts (e.g. 
Springel 2005 , Tasker et al. 2008 ). Tests of the formation of cosmo- 
logical haloes assuming non-radiative hydrodynamics have demon- 
strated that SPH produces similar results to those of Eulerian adap- 
tive mesh refinement (AMR) schemes beyond the central 10 per- 
cent of the halo virial radius ( |Frenk et al.j 1999, Kravtsov et al. 
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|2005||Voit et al.|2 005 ). This agreement in the limit of non-radiative 
physics therefore suggests that, in spite of its inability to resolve 
shocks well, SPH schemes model gravitational shock heating ade- 
quately. 



2.2 Halo and galaxy identification 

As decribed in C09, we identify bound haloes using the SUBFIND 
algorithm presented by Dolag et al. (2009), which extends the stan- 
dard implementation (Springel et al. 2001) by also considering 
baryonic particles when identifying self-bound substructures. This 
procedure first finds dark matter haloes using a friends-of-friends 
(FoF) algorithm, with the standard linking length in units of the in- 
terparticle separation (b = 0.2; Davis et al. 1985). It also associates 
all baryonic particles with their nearest neighbour dark matter parti- 
cle. The aggregated properties of baryonic particles associated with 
grouped dark matter particles define the baryonic properties of each 
FoF group. Halo substructures are then identified using a topologi- 
cal unbinding algorithm. This provides an unambiguous definition 
of a galaxy within the simulations, namely the set of star particles 
bound to individual subhaloes. The gas bound to each subhalo then 
forms the interstellar medium (ISM) and the hot corona. Since in- 
dividual haloes may have more than one 'subhalo', they may host 
more than one galaxy. In analogy to semi-analytic models, the stars 
associated with the most massive subhalo of a FoF group are hence 
defined as central galaxies, whilst stars associated with substruc- 
tures are satellites. 



2.3 Morphological classification 

The resolution of our simulations is sufficient to allow a mor- 
phological classification of the galaxies into disc- and spheroid- 
dominated types, based on their dynamics. We assume a simple 
two-component model: i) a dispersion-supported spheroid, and ii) a 
rotationally- supported disc. The centre of (baryonic + dark) mass of 
the main subhalo of each FoF group is used as the starting point for 
an iterative procedure that computes the centre of mass of all star 
particles within a sphere that, at each iteration, is decreased in ra- 
dius by 2 per cent and re-centred on the new centre of stellar mass. 
The procedure concludes when the sphere encloses fewer than 32 
star particles; we have checked that this scheme yields robust galac- 
tic centres. The smaller of three times the stellar half mass radius, 
r\/2, or the distance to the furthest bound star particle is then used 
as a 'boundary' for the galaxy, r ga i. This is to exclude the contri- 
bution of the diffuse 'intrahalo stars' that SUBFIND associates with 
the potential of the most massive halo substructure. This is an im- 
portant practice, since our sample includes a number of massive 
galaxies for which up to 30 per cent of the bound stellar mass is 
(by this definition) in the form of intrahalo stars. 

To decompose the stellar mass of galaxies into spheroid and 
disc components, we use a procedure based on the method intro- 
duced by Abadi et al. (2003). We compute the angular momentum, 
L, of the N t star particles within r ga i: 






I"i X Pi 



(1) 



where r is the radial vector with respect to the centre of the system 
and p is the linear momentum vector corrected for the bulk pecu- 
liar velocity of the system. The assumption that the spheroid is fully 
dispersion supported requires that it should have no net angular mo- 
mentum, and so the spheroid mass can be reasonably approximated 
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Figure 1. Histogram of disc-to-total stellar mass ratios (D/T) for the five 
GIMIC regions. Results from the high-resolution — 2a region are also shown 
to illustrate the sensitivity of morphology to resolution. Approximately half 
of all galaxies in the simulations are disc-dominated (i.e. D/T > 0.5) at 
z = 0. To create a sample that mimics "late-type" galaxies we adopt a 
criterion of D/T > 0.3 for inclusion in our sample. 



as twice the summed mass of particles that counter-rotate with re- 
spect to L. The remaining stellar mass of the galaxy then comprises 
the disc. 

The distribution of D/T for all galaxies with M* > 10 10 M Q 
is shown, for the five GIMIC regions, in Fig. [T] The total number 
of galaxies included from the five intermediate resolution regions 
is 1267, whilst 111 galaxies are included from the high-resolution 
—2a region. In each region, approximately half of all galaxies are 
disc-dominated (D/T > 0.5); this result is consistent with the mor- 

■ Driver 



phological analyses of the Millennium Galaxy Catalogue by 



et al. (2006) and of the SDSS catalogue by |Benson et al.| l 



2007). 



Since we wish to create a sample that excludes only obviously "el- 
liptical" galaxies, we adopt a slightly lower threshold and define 
"disc galaxies" as those with a disc-to-total stellar mass fraction 
> 0.3. The precise choice of this fraction is unimportant for the 
purposes of this study; we have checked that our results are not 
affected by adopting a value of 0.2 or 0.5 instead. 

The high proportion of disc-dominated galaxies is a notewor- 
thy success of our simulations. It is a well-known problem of hy- 
drodynamical simulations of this type that excessive cooling in 
small halos, allied to angular momentum transport during merg- 
ers, results in discs with much lower angular momentum than ob- 
served for real galaxies (Weil et al. 1998). Recently, high-resolution 
simulations of individual objects have demonstrated the importance 
of feedback in mitigating this problem (e.g. Okamot o et al.|2005} 
Scannapiec o et al.|2008fr , a finding supported by the large sample 
of galaxies we present here. We plan to investigate this interesting 
result in future work. 



2.4 Computing luminosities 

Optical and X-ray luminosities are computed in postprocessing. 
The former are calculated individually for star particles, consider- 
ing them as simple stellar populations (SSP). The masses of in- 
dividual stellar particles produced in our simulations (~ 10 6 — 
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Figure 2. Synthetic surface brightness map at 2 = of a galaxy drawn 
from our sample of GIMIC galaxies. The map is 175 kpc in each dimension, 
projected over the same depth. The galaxy is edge-on to the projection axis. 
The colour mapping shows the surface brightness of soft X-ray emission, 
whilst the blue contours trace the underlying i\"-band surface brightness of 
the galaxy. 



10 7 Mq) are comparable to, or greater than, the typical mass of star 
clusters, so it is reasonable to assume that they can be described as 
an SSP. We assume a universal initial mass function (IMF; C habrier| 
2003) and store the age and metallicity of star particles. Thus a 
spectral energy distribution (SED) for each particle can be derived 
by interpolation over the GALAXEV models of Bruzual & Chariot 
(2003). The optical luminosity of each star particle is obtained by 
integrating the product of its SED with the appropriate filter trans- 
mission function; the overall broadband luminosity of a galaxy is 
then defined as the sum of the luminosities of all star particles 
within r ga i. Since we are primarily concerned here with AT-band 
luminosities, we do not expect dust extinction to be important. 

Gas phase X-ray luminosities are also computed on a per- 
particle basis. We arrive at the overall X-ray luminosity of a galaxy 
by summing the luminosities of all gas particles bound to its sub- 
halo. The X-ray luminosity of the j gas particle is computed as: 



j x,j 



% e ,jnn,jAjVj 



(2) 
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[X C (Z 3 ) + XiiZ^^Z^mn KZj)m H J > 

where p is the particle gas density, m gaa is the particle mass, the 
volume is V = m gas /p, n c , «h and m are the number densities 
of electrons, hydrogen, and ions, respectively, X c = n c /nu, Xi = 
rii/rin, Z is the metallicity, p is the mean molecular weight, ttih is 
the mass of a hydrogen atom, and A is the cooling function in units 
of ergs cm 3 s _1 (integrated over some appropriate passband, such 
as 0.5-2.0 keV). We compute A by interpolating a pre-computed 
table generated using the Astrophysical Plasma Emission CodJj 



(APEC, v 1.3.1, see |Smith et al.|200l) under the assumption that the 
gas is an optically thin plasma in collisional ionisation equilibrium. 
APEC cooling rates are computed on an element-by-element basis 
and summed to yield the total cooling rate of each particle, i.e., 



Ai(T J -)=£)A i ,*(T i ; 



(3) 



where T is the gas temperature and \j,k(Tj) is the cooling func- 
tion for element species k for the j th particle. The summation is 
performed over the 1 1 most important elements for cooling (H, He, 
C, Ca, N, O, Ne, Mg, S, Si, Fe), which are individually and self- 
consistently tracked during the simulation. APEC assumes the solar 
abundance ratios of Anders & Grevesse (1989), but it is straight- 
forward to modify the spectra for arbitrary abundances, which we 
have done. Note that our scheme automatically excludes gas within 
the cool disc of the interstellar medium (ISM, i.e. those for which 
nu > 0.1 cm' 3 ) since we assign these particles a temperature of 
T = 10 4 K, which is below the minimum temperature for which 
APEC returns a non-zero X-ray luminosity. Therefore we consider, 
by construction, only extra-planar emission. 

It should be noted that the APEC cooling tables assume pure 
collisional ionisation equilibrium and hence neglect the extragalac- 
tic UV/X-ray background. However, the effect of photoheating on 
the derived X-ray properties is small in the regime we explore here. 
Wiersma et al. (2009a) show that the difference in cooling time for 
hot (T > 2 x 10 5 K) plasmas under collisional and photoioni- 
sational equilibrium, with and without photoheating, is essentially 
negligible at all densities in the absence of metals, but for low- 
density plasmas (bh < 10~ 4 cm -3 ) enriched to solar abundance, 
the cooling time can be a factor of > 2 greater in the presence 
of a metagalactic photoionising UV/X-ray flux. However, we find 
that only a small fraction of the total X-ray luminosity of L* haloes 
comes from gas at such low densities. It is important to note though 
that this does not imply that photoheating is unimportant in the for- 
mation and evolution of disc galaxies. Indeed, it is possible that the 
progenitors of these systems could have been substantially affected 
by the UV/X-ray background. 



2.5 The sample of simulated galaxies 

In the following section, we compare the results of the GIMIC sim- 
ulations with observational measurements. The galaxies are drawn 
from all 5 intermediate-resolution GIMIC regions, and we impose 
a number of selection criteria in order to construct a sample simi- 
lar to observational samples. We select central (i.e. the most mas- 
sive galaxy within a FoF halo), disc galaxies (D/T > 0.3) whose 
stellar mass lies in the range 10 10 < M t < 10 11 ' 7 M Q , roughly 
corresponding to L* at z — 0. Additionally, we require that the 
subhalo with which each galaxy is associated should account for 
at least 90 per cent of the mass of the parent FoF group, i.e. 
Afsub > 0.9A/f o f. This criterion is intended to select only isolated 
systems by excluding galaxies that are interacting or are members 
of galaxy groups and clusters. The 5 regions yield a sample of 458 
galaxies fulfilling these criteria. 

To illustrate the correspondence of our methods with obser- 
vational techniques, we show in Fig. [2] a synthetic surface bright- 
ness map of the diffuse soft X-ray and A'-band optical emission, 



1 To maintain strict consistency with the implementation of radiative cool- 
ing in the simulations, it would be more appropriate to use cooling rates pre- 
dicted by the CLOUDY software package (Ferland et al. 1998). The APEC, 



however, is more widely used in the analysis of X-ray data, which is why we 
have adopted it here. We verified that using CLOUDY instead gives nearly 
identical results. 
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Figure 3. The present-day 0.5-2.0 keV X-ray luminosity as a function of K- 
band luminosity. Simulated galaxies are represented by black dots. For the 
observational data, extra-planar point source-corrected X-ray luminosities 
have been extracted from a number of recent studies (see legend and text). 
The corresponding A'-band luminosities of these galaxies were extracted 
from the online 2MASS database at IPAC. The simulations reproduce both 
the scaling and the scatter of the observations. 



computed for a galaxy drawn from our sample. The galaxy has 
the following properties: M200 = 1O 12O M , M* = 1O 1O6 M , 
L K = 10 10 - 8 L K ,©,Lx,o.5-2.okcV = 10 391 erg s" 1 . In general the 
surface brightness increases towards the galactic centre, in keeping 
with the WF91 picture, but close to the galaxy disc the X-ray emis- 
sion structure becomes more complex, highlighting the necessity 
of employing hydrodynamical simulations to calculate the detailed 
evolution of the gas. Extra-planar emission is clearly visible, with 
the contours tracing the dense disc of the ISM at small radii and be- 
coming more circular (in projection) at larger radii, except where 
their morphology is disturbed by accreting substructure. 



3 X-RAY CORONAE OF DISC GALAXIES: 
COMPARISON WITH OBSERVATIONS 

In this section we examine the X-ray scaling relations found in our 
hydrodynamic simulations. In Fig. [3] we plot the soft (0.5-2.0 keV) 
X-ray luminosity of the simulated galaxies (black dots), as a func- 
tion of their i\"-band luminosity. The colour symbols are measure- 
ments for real disc galaxies collated from the literature {Benson 
eTaTp OOO; Stricklan d~et~aTl|2004] [Wa ng 200 5] |Tullmann et al 
20061 |Li et al.||2007| |Owen & Warwick||2009| |Rasmussen eTal 



2009|[Sun et al.|2009) . Systems for which X-ray detections were 
made are denoted by solid symbols, whilst those with only upper 
limits are denoted by open symbols with down arrows. In the case 
of X-ray detections, we use only the reported extra-planar X-ray lu- 
minosity, which has been corrected for point source contamination 
(by spatial filtering of bright sources and the inclusion of a power- 
law component in the spectral fitting for faint sources). Where X- 
ray luminosities have been quoted in other passbands, we correct 
into the 0.5 — 2.0 keV band using the PlMMS softwar^](Mukai 



1993). K-banA luminosities have been extracted from iPAC's on- 
line database for the Two Micron All-Sky Survey (2MASS, [Skrut-| 
|skie et al.|20Qg) . 

The K-b&nA luminosity is a good tracer of stellar mass, which 
itself is expected to correlate with halo mass. The simulations yield 
an X-ray vs. if-band luminosity relation that is in very good agree- 
ment with the observed detections, and is consistent with the upper 
limits for non-detections. It is interesting that the simulations also 
reproduce the large scatter in the observations; we discuss the ori- 
gin of this scatter in § |3.1| It should be noted that we have not tuned 
any element of the simulation subgrid physics, for example, the 
supernovae wind velocity or mass-loading, in order to reproduce 
the observed relation. As discussed in C09 (see also Schaye et al. 
2010), the adopted feedback parameters were chosen only so as to 
produce a reasonable match to the peak of the overall star formation 
history of the Universe at z ~ 2 — 3. 

The canonical picture posits that the soft X-ray luminosity of 
hot coronae increases steeply with the overall mass of the system 
(i.e. the galaxy and its host dark matter halo). It should therefore 
be easiest to detect soft X-ray emission from the most massive sys- 
tems. The virial mass of a halo is commonly approximated by the 
mass, M200, of the sphere of radius r2oo about the centre of the 
potential that encloses a mean density of 200 times the critical den- 
sity for closure. In practice, estimating the total mass is difficult 
and requires recourse to observational proxies. For disc galaxies, 
the rotation speed, v rot , is often used since, under the assumption 
that the halo has an isothermal density profile, p oc r -2 , the cir- 
cular velocity profile, v c (r) — [GM (< r)/r] 1 ' 2 , is simply a con- 
stant at all radii. This enables the rotation speed of a centrifugally 
supported disc to be equated to the virial circular velocity of the 
halo, «2oo = ^c(»"2oo)- However, if the structure of the galaxy-halo 
system deviates from the isothermal profile, the circular velocity 
profile is not constant, and the relationship between the rotation 
velocity of a galaxy and the virial mass of the galaxy's host halo 
becomes more complicated. 

The Navarro-Frenk-White (hereafter NFW, Nav arro et al.| 
|1995| |1996| |1997| > profile, shown to be a near universal density 
law for dark matter haloes, deviates from the isothermal form at 
both small and large radii (where asymptotically p oc r _1 and 
p oc r -3 respectively). Thus, the circular velocity curve v c (r) 
of the NFW profile varies with radius. For example, an NFW 
halo with concentration c = 10 exhibits a peak circular velocity, 
Umax = max[u c (r)] = 1.2n2oo at r ~ 0.2r 2 oo- The density pro- 
file of galaxy haloes will be modified by the baryons cooling in 
the halo (e. g. |Blumenthal et al.|1986||Jesseit et al.|2002"l |Gnedin| 



http://heasarc.nasa.gov/docs/software/tools/pimms.html 



|et al.|2004|[Abadi et al.|2009[|Duffy et al.|2010) . In general, cold 
baryons will be preferentially deposited at the halo centre leading 
to a steepening of the central gravitational potential (both because 
of the presence of the baryons and the adiabatic contraction of the 
dark matter halo), raising the peak of the circular velocity curve to 
«max/«200 ~ (1.5 — 2.0). These simple theoretical arguments in- 
dicate that the halo mass cannot be accurately estimated by equat- 
ing the rotation speed of galaxy discs with l>200- Both modelling 
improvements - switching from isothermal to NFW profiles, and 
accounting for the effects of baryons - are important, with the latter 
perhaps being somewhat more so. 

In Fig.Elwe plot the relationship between the soft X-ray lumi- 
nosity and two different 'measures' of the disc rotation speed for 
the simulated galaxies. In the left-hand panel, we identify v TO t with 
the circular velocity at the virial radius (as in the isothermal pro- 
file); in the right-hand panel we use instead w max computed from 
all mass components, i.e. gas, stars and dark matter. The obser- 
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Figure 4. The present-day 0.5-2.0 keV X-ray luminosity as a function of disc rotation velocity. In the left-hand panel, we identify the disc rotation speed with 
"200 (the halo circular velocity at r2oo). following WF91's assumption of an isothermal halo density profile. In the right-hand panel, we use instead " max , the 
maximum of the halo radial circular velocity profile. The observational data (various coloured symbols) are the same in both panels; the disc rotation velocities 
were extracted from the HYPERLEDA database. Black dots represent galaxies drawn from our simulations. We give the associated mass, M200> f° r a given 
"200 on the upper axis of the left-hand panel, since the conversion between the two is trivial. The data show that the assumption v TO t = "200 is inaccurate; 
the approximation v Iot = " max produces much better agreement between the simulations and the observational data. 



vational data are the same in both panels; rotation velocities are 
inclination-corrected 21 -cm measuremen ts extracted from the Hy- 
PERLeda databasfl JPaturel et al.|2003 i. 

When we adopt the inaccurate approximation, v Iot ~ U200, 
the simulated galaxies are much too X-ray luminous at a given 
value of w rot compared to the observational data. A similar con- 
clusion was drawn by BOO and others when comparing the analytic 
predictions of WF91 to observational data. However, when we use 
Umax as a measure of u ro t the simulations are in much closer agree- 
ment with the data. This is a much more realistic assumption since 
r m ax is typically < lOkpc, comparable to the sizes of the simu- 
lated discs*! The right-hand panel of Fig. kl immediately indicates 
that a significant improvement can be made to analytic estimates 
of Lx simply by adopting a realistic density profile for the halo. 
However, this effect alone is insufficient to reconcile the analytic 
predictions of WF91 with the observations and, as we discuss be- 
low in §4, a reduction in the density of the hot gas halo is also 
required. 

In several cases where diffuse X-ray emission has been de- 
tected from disc galaxies, a correlation between the X-ray lumi- 
nosity of the corona and the star formation rate of the galaxy has 
been noted. Furthermore, there exist a few spectacular cases, such 
as the starburst galaxy M82 ( |Strickland et al.|2004] |Strickland &| 
Heckman 2007), where the detected X-ray emission has a bicon- 



3 http://leda.univ-lyonl.fr/ 

4 In principle, the rotation speed of the discs could be calculated directly 
in the simulations. In practice, this is inaccurate because of small particle 
number and resolution effects. We defer a discussion of the structure of the 
discs to a more detailed paper. 
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Figure 5. The present-day 0.5-2.0 keV X-ray luminosity as a function of 
star formation rate. The star formation rates of observed galaxies have been 
derived from near-infrared luminosities, using the relation of Kennicutt 
1 19981, corrected to our assumed Chabrier 1 2003 1 IMF. Note that this is 
a small correction (see text). The simulated galaxies, shown as black dots, 
are consistent with both the scaling and scatter of the data. 



ical morphology. On this basis, several authors have concluded 
that the X-ray emission must originate from hot outflowing gas 
that has recently been heated by supernovae, rather than from the 
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gravitationally-heated inflowing gas of a cooling flow. In Fig. B] 
we show the relation between the soft X-ray luminosity and the 
star formation rate (SFR), M*, of our simulated galaxies. Over- 
plotted are estimates of AL for several observed systems, obtained 
from far infrared (FIR) continuum luminosities using the relation 
of Kennicutt ( 1998). For consistency, we correct these SFRs to the 
Chabrier ( 2003 1 IMF assumed in the simulations rather than that 
adopted by Kennicutt. For this, we assume that our galaxies ex- 
perience ongoing star formation, and hence that the FIR continuum 
luminosity is dominated by stars more massive than 5 Mq. We then 
scale the inferred SFR by the mass fraction of each IMF comprised 
by these massive stars (0.23 & 0.29 for Kennicutt and Chabrier, re- 
spectively). For the same FIR continuum luminosity, the inferred 
SFR corresponding to the Chabrier IMF is therefore 79 per cent of 
that for the Kennicutt IMF, a relatively small correction. 

Whilst a correlation between X-ray luminosity and star for- 
mation rate has been found in previous studies that focussed ex- 
clusively on galaxies with high X-ray luminosities (e.g. Strickland 
|et al.|2004| >, it is interesting to note that the addition of measure- 
ments for low X-ray luminosity systems in Fig.[5]renders the con- 
nection much less compelling, even if one excludes outlying sys- 
tems such as Ml 04. The simulations are broadly compatible with 
these observations (for M* < 5 Mq yr _1 ), exhibiting a loose cor- 
relation between Lx and AL* with a scatter in Lx of more than a 
factor of 10 at fixed A/*. Interestingly, the simulations indicate an 
additional reason for caution when interpreting the perceived cor- 
relation between Lx and M*, in that both quantities are expected 
to scale with the overall mass of the system. Any correlation be- 
tween these quantities does not, by itself, demonstrate that the X- 
ray emission originates in outflows driven by star formation, and 
we challenge the generality of this frequently made interpretation 
in § [4311 

3.1 Scatter in the Lx Lk relation 

Finally, we consider the system-to-system scatter. The observa- 
tional scalings show considerable scatter, for example in the Lx — 
Lk, Lx — v rot and Lx — M* relations, which is also present in 
the simulations. To investigate the origin of this scatter, we plot the 
scaling between X-ray luminosity and it -band luminosity, stellar 
mass and virial mass, in Fig. [6] As anticipated, the Lf-band lu- 
minosity traces the stellar mass of the galaxy very closely. These 
panels demonstrate that the scatter in the X-ray scaling relations 
stems primarily from the scatter in stellar mass at fixed halo virial 
mass, which itself is sensitive to the assembly history of the galaxy. 
By contrast, the scatter in Lx at fixed halo mass is relatively small 
for all halo masses considered here, raising the interesting possibil- 
ity that X-ray properties might be understood in greater detail by 
focusing on systems for which a more accurate and precise mea- 
surement of the halo mass can be found, such as those derived 
from strong lensing l |Ferreras et al.|2005| >, planetary nebulae (Ro- 
|manowsky et al.| 2003 ; Napolitano et al. 2009) or the properties of 
globular clusters ( Spitler & Forbes 2009 ). 

The scatter in the Lx — Lk plane for elliptical galaxies has 
received considerable attention, since much larger samples of el- 
lipticals can be compiled. The standard interpretation of the scatter 
in Lx at fixed Lk in these samples is, as seen for our sample of 
disc galaxies, that such galaxies are hosted by dark matter haloes 
drawn from a wide range of total masses (i.e. AL200X leading to a 
similarly wide range of X-ray luminosities (Mathews & Brighenti 
|2003| [Mathews et al. 2006). More recently, an analysis of archival 
Chandra data by Jeltema et al. (2008) uncovered a significantly 



greater fraction of X-ray detected Lk > L* galaxies in groups 
relative to clusters. They concluded that environmental effects, for 
example ram pressure stripping, are a considerable source of scatter 
at the bright end of the Lx — Lk relation. Since we have focussed 
here on isolated L* disc galaxies, the coronae in our sample will 
not have been affected by ram pressure stripping. Nonetheless, we 
find considerable scatter in the Lx — Lk plane, indicating that some 
fraction of this scatter is due to the range of M* at fixed A/200 (and 
vice- versa). This is illustrated succintly by the contrasting scatter 
in the Lx — AL and Lx — A/200 relations (Fig. [6] centre and right- 
hand panels respectively). We note, however, that this does not pre- 
clude an effect of the kind found by Jeltema et al. (2008 ) for the 
early-type galaxies hosted by massive clusters in these simulations. 



4 THEORETICAL INTERPRETATION 

Having established that the GIMIC simulation generally reproduces 
the basic scaling relations observed for the soft X-ray luminosity 
of disc galaxies, it seems unlikely that a fundamental shortcoming 
in the canonical theory of disc galaxy formation is the reason for 
its large overprediction of the X-ray luminosities of galaxies. In- 
stead, it seems more likely that inaccuracies in one or more of the 
approximations adopted by WF91 are to blame. We explore these 
approximations in this section. 



4.1 The analytic basis for X-ray coronae 

In sufficiently massive haloes it is expected that accreting gas will 
be shock-heated to a temperature that reflects the depth of the dark 
matter potential well, i.e. the halo virial temperature, 



T 20 o = \ ^4,0 = 3.6xl0 5 K 
3 Kb 



M \ / «200 

/ V: 



0.59 



lOOkms- 



•(4) 



In massive haloes, the post-shock cooling time of the gas will, in 
general, be greater than the dynamical time of the halo, leading 
to the formation of a quasi-hydrostatic corona. Cooling of some 
fraction of this hot gas leads to a cooling flow from which a ro- 
tationally supported gas disc forms. A stellar disc then grows from 
quiescent star formation within the cold gas disc. This scenario was 
first put forward by White & Rees ( 1978 ) and was implemented in 
the CDM cosmology by WF91, leading to the important conclu- 
sion that massive galaxies should, at the present day, be embed- 
ded within extended atmosheres of hot, X-ray luminous gas. Disc 
galaxies represent the cleanest observational test of this picture be- 
cause, according to the WF91 model, they are systems with a rel- 
atively quiescent star formation history for which the central as- 
sumptions of the model appear most justified. These assumptions 
are as follows: 

• WF91 assumed that both the dark matter halo and the gas 
corona follow an isothermal density profile, p oc r~ 2 . The gas 
profile is normalised to obtain a desired overall gas mass fraction 
within the virial radius. In the simplest case, where cooling and star 
formation are inefficient and feedback is energetically incapable of 
ejecting gas from the halo, / gas (f2oo) ~ fib/f2 m . Alternative an- 
alytic forms (such as the NFW profile) have also been explored in 
the context of the WF91 formalism (see, e.g., BOO). 

• The temperature of the gas is equal to the halo virial tempera- 
ture, T 20 o- 

• In the simplest case, the gas is assumed to have a primordial 
composition of hydrogen and helium (Y = 0.24). The inclusion of 
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Figure 6. The scaling of soft X-ray luminosity with i\~-band luminosity (left), stellar mass (centre) and halo mass (right). The similarity of the left-hand and 
centre panels demonstrates that the ftT-band luminosity (computed using the GALAXEV model of Bruzual & Chariot 2003 1 traces the stellar mass very closely, 
with a nearly 1 : 1 relation between Lk ™d M* ■ The contrast in the scatter exhibited in the centre and right-hand panels indicates a significant scatter in stellar 
mass at fixed virial mass (consistent with the top right panel of Fig. 17}, and it is this effect that drives the scatter in the Lx — Lk relation. 



radiative emission due to metal lines increases the predicted X-ray 
luminosity. 

• The gas cools in a single phase quasi-hydrostatic cooling flow 
at a rate, M coo i, that is determined by the assumed gas density, tem- 
perature, and metallicity distributions (see Eqns. 20-22 of WF91). 
This leads to an X-ray luminosity, Lx oc Mcooi^oo (as first high- 
lighted by |Thomas et al.|1986) . 

These assumptions enable the complex physical processes un- 
derlying galaxy formation to be condensed into a relatively simple 
set of analytic equations. We expect that, to some degree, all of 
these assumptions will be violated in the real world and in cos- 
mological hydrodynamic simulations. In the following subsection, 
we quantify to extent to which the analytical approximations are 
violated in the simulations and the impact on the predicted X-ray 
emission from massive disc galaxies. 



4.2 Why are the predicted X-ray luminosities now in 
agreement with observations? 



Toft et al. (2002) concluded that the X-ray luminosity of galaxies in 
their hydrodynamic simulations was lower than in analytic models 
because efficient radiative cooling over cosmic history substantially 
reduced the hot gas fraction of dark matter haloes by z = 0. On the 
other hand, C09 demonstrated that supernova winds strongly re- 
duce the mass of dense, highly X-ray luminous gas near the centres 
of haloes with circular velocity ^200 J; 300 km s _1 at z — 0. 
Thus, both cooling and heating processes can act to reduce the X- 
ray luminosity of galaxies. 

We begin by scrutinising the halo baryon content of our galaxy 
sample in greater detail, in Fig. [7] The first three panels show the 
overall baryon fraction (gas + stars, top left), the stellar fraction (top 
right) and the gas fraction (bottom left), each normalised by the cos- 
mic baryon fraction. The final panel (bottom right) shows the frac- 
tion of halo gas that is X-ray luminous (i.e. T > 2 x 10 5 K). As 
a consequence of supernova feedback, the overall baryon fraction 
is a strong function of halo mass: haloes with M200 ;$ 10 12 Mq 
have overall baryon fractions that are typically half of the cos- 
mic value. At higher masses, halo baryon fractions approach the 



/b — 0.9f2b/f2 m limit found by Frenk et al. j!996f and Crain et al. 
(2007) in the non-radiative regime, indicated in the top left panel 
by a dotted line. The most massive haloes in our sample, whose 
masses approach those of small galaxy groups, are therefore es- 
sentially 'baryonically closed' systems^] However, for all systems 
(including the most massive ones) a non-negligible fraction of the 
baryons that reside within the virial radius are in the form of stars 
(top-right panel of the figure) and, in the case of low mass haloes, 
cold gas (lower panels). The stellar mass fraction in haloes with 
mass M200 < 10 12 Mo is, on average, similar to the fraction found 
by |Guo et aL] {2010| l for galaxies in the Sloan survey; for larger 
haloes, the stellar masses are about a factor of two too large, but on 
these scales our neglect of AGN feedback could be significant. 

The low mass haloes in our simulated sample have baryon 
mass fractions that have been reduced relative to the cosmic mean 
by feedback. The mean gas density in these haloes is also lower 
than assumed in the canonical picture. This deficit of hot gas im- 
mediately points to a likely cause for the lower X-ray luminosities 
found here compared to WF9 1 . However, even those systems that 
we identified as baryonically closed, and that therefore still retain 
relatively high mean densities, exhibit significantly lower X-ray lu- 
minosities than predicted analytically. Such low X-ray luminosities 
warrant an explanation. 

A clearer understanding can be obtained from radial profiles 
of the thermodynamic state of the hot gas, rather than just the mass 
fraction (since X-ray luminosity exhibits a strong dependence on 
the density distribution and is particularly sensitive to its profile 
at the halo centre). Spherically-averaged radial profiles of the key 
thermodynamic quantities at z — for all gas with T^2x 10 5 K 
are shown in Fig. [8] density (top left), temperature (top right), en- 
tropy (bottom left) and cooling time (bottom right). The profiles are 
split into halo mass bins and the plotted lines represent the median 



5 We note here that this result violates observational constraints on the 
galaxy group scale derived from a combination of optical and X-ray mea- 
sures, but this can be understood in terms of the much greater role expected 
to be played by feedback from active galactic nuclei (AGN) in this high- 
mass regime (e.g. McCarthy et al. 2009 1. We do not attempt to model this 
particular feedback mechanism in these simulations. 
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Figure 7. The baryon content, at z = 0, of the haloes included in the simulated galaxy sample. The total baryon fraction, /^ = (M* + M gas )/Af200 (top 
left), the stellar fraction /* = M+ /M20Q (top right), and the gas fraction / gas = M gas /A/200 (bottom left) are normalised by the cosmic baryon fraction 
of fl^/flm. The bottom right panel shows the mass fraction of gas that is hot (T > 2 X 10 5 K), and thus potentially X-ray luminous. The dotted line in the 
top-left panel marks the baryon fraction found by Frenk et al. 1 1996i and Crain et al. 12007 1 in the case of non-radiative gas. 



profiles of haloes in each bin. The profiles are discussed in turn 
below. 

• Density. The gas density profiles of the simulated haloes devi- 
ate significantly from the favoured analytic forms, isothermal (dot- 
ted line) and NFW (dashed line). The least massive galaxies ex- 
hibit the greatest deviation from the analytic models. Interestingly, 
at large radii where the cooling time of the gas is relatively long 
(see bottom right panel), the profiles show a monotonic trend in 
mass, reflecting an increasing deviation from the analytic profiles in 
less massive haloes. Suppression of the hot gas density, by cooling 
(or more specifically, star formation, since this is the terminal state 
for cooled baryons) and feedback, is therefore most pronounced in 
these systems. Most importantly, the profiles indicate that star for- 
mation and feedback preferentially suppress the hot gas density at 
the centre of haloes, and it is the central density to which the X-ray 
luminosity is most sensitive. We note that, typically, half of the to- 
tal X-ray luminosity of a simulated galaxy comes from gas within 
the central 10 per cent of the virial radius. 

• Temperature. The gas at the halo centre is not isothermal, 



but shows a mild negative temperature gradient. Gas temperatures 
reach up to 3 times the virial value in the central regions and fall 
slightly below it at the virial radius. The origin of central super- 
virial temperatures is twofold: firstly, heating by feedback and, sec- 
ondly, gravitational compression that results from the steepening of 
the halo's central potential as it accumulates cold baryons. While 
this illustrates that the assumption that the gas is at the virial tem- 
perature is incorrect in detail, the error incurred by adopting this 
approximation is not large. Typically, we find that it affects the lu- 
minosity by no more than 30 per cent, and it exhibits a mild de- 
pendence on the heavy element abundances. This stems from the 
fact that the virial temperature of a halo with M200 ~ 10 12 Mq is 
~ 10 6 K, which roughly corresponds to the minimum of the cool- 
ing function. Thus, a change in T by a factor of a few yields only a 
small change in A(T, Z). 

• Entropy. The density and temperature profiles reflect the un- 
derlying entropy configuration and the depth of the gravitational 
potential well (which is dominated by dark matter). Unlike the den- 
sity or temperature, which are sensitive to processes that compress 
or rarify the gas, entropy is conserved in any adiabatic process. 
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Figure 8. Present-day luminosity weighted, three-dimensional, spherically averaged radial profiles of the hot gas bound to the haloes of disc galaxies in our 
sample. Galaxies are binned in log 10 M200 at 2 = (see legend, units are M©). Shown are the density (top left), temperature (top right), entropy (bottom 
left), and cooling time (bottom right). The density profiles clearly deviate from the isothermal (dotted line) and NFW (dashed line) profiles, commonly adopted 
in analytic models. Feedback from supernova-driven winds drives gas from the dense halo centre, strongly suppressing the X-ray luminosity, particularly for 
less massive systems. Feedback and non-gravitational heating also render the temperature profile super-virial. The net effect is that the entropy profiles deviate 
significantly from the self-similar scaling expected in the non-radiative regime, shown as a dotted line in the bottom-left panel. This entropy elevation lengthens 
the cooling time relative to that expected from analytic modelling under the assumption of primordial composition gas. This is seen in the bottom-right panel 
where the dotted lines correspond to the canonical WF91 model as calculated by Benson et ah] (2000). 



Non-adiabatic cooling and heating processes are therefore most 
clearly described by entropy, since cooling always reduces the local 
entropy of a system, whilst heating always raises it. 

Entropy therefore maintains a record of the thermodynamic his- 
tory of the gas (e.g. |Voit et al.|2003|[Voit|2005| l. We consider here 
the'virial entropy', which is denned as 



£200 



feTioo 



(5) 



K,200] 2 / 3 ' 

As can be clearly seen in the bottom-left panel of Fig. [5] the hot 
gas in our simulated haloes exhibit s a very strong de viation from 
the self-similar entropy scalinjfM of Voit et al. 1 2005 [ dotted line) 



which is the entropy distribution expected in the case where cool- 
ing and feedback are negligible. The 'excess entropy' with respect 
to the self-similar scaling - for which there is already firm observa- 
tional evidence on the scale of systems with M200 > 10 Mq (e.g. 
|Balogh et al.[T999l| Johnson et al.|2009[|Sun et al.|2009) - therefore 
represents a strong indication that non-gravitational processes are 
vitally important in establishing the thermodynamic properties of 
hot galactic coronae. 

• Cooling time. As a result of the increased entropy (and, in turn, 



The dotted line represents a power-law fit to the entropy profiles of a 



sample of galaxy groups and clusters simulated by |Voit et al.H2005) with 
non-radiative hydrodynamics. Over the range 0.1 — 1.0r2oo> hydrodynamic 
simulations employing either adaptive mesh refinement (AMR) or SPH al- 
gorithms yield consistent profiles. 
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Figure 9. Present-day three-dimensional, spherically-averaged cumulative 
radial profile of all mass (i.e. gas, stars and dark matter) comprising the 
main subhalo of disc galaxies included in our sample. As per Fig. [8] galax- 
ies are binned in log 10 M200 at 2 = 0. The profiles exhibit excellent agree- 
ment with the NFW (dashed line) profile at large radii, but are more con- 
centrated than both the isothermal (dotted line) and NFW profiles within 
~ 0.3r200 due to the accumulation of cooled baryons. 



reduced gas density), the cooling time of the hot gas in the simu- 
lated galaxies is longer than in the WF91 model. Therefore, if cool- 
ing dominates non-gravitational heating and there is a net inflow of 
hot gas (see §4.3 below), the relatively long cooling time of the gas 
compared to the halo dynamical time (typically ~ 1 — 2 Gyr) im- 
plies that the flow is quasi-hydrostatic, as assumed by WF91. We 
note that the cooling times we compute for the simulated galax- 
ies are self-consistent, in that they account for the thermodynam- 
ics and chemical abundances of gas particles. (Typically, we find 
the mean emission-weighted metallicity of the hot gas to range be- 
tween 0.3 — 1.3 Zq.) The profiles are compared with the analytic 
counterparts proposed by WF91 for isothermal gas of primordial 
composition (dotted lines). The discrepancy between the simula- 
tions and the WF91 model would be greater if metal line emission 
was taken into account in the latter. 

For reference, we show in Fig. [9] the cumulative total mass 
radial profiles (i.e. including gas, stars and dark matter) of the 
selected systems, again as a function of halo mass. The profiles 
trace the NFW form (dashed line) closely at large radii, but for 
r < 0.3r2oo they are more concentrated than both the NFW form 
(for a reasonable choice of the concentration parameter) and the 
isothermal form (dotted line), down to the innermost radius reli- 
ably probed given the resolution of the simulations, which can be 
approximated as ~ 4 smoothing lengths (Hernquist & Katz 1989), 
corresponding to a few percent of r2oo- The boost in concentration 
for the central mass profile follows from the accumulation of cooled 
baryons that also gives rise to the mildly super-virial gas tempera- 
tures seen in Fig. [8] and is broadly consistent with the picture of 
baryon-induced halo profile backreactions recently discussed by 
|Duffy et al.lpOTOl l. 

To summarise, the X-ray luminosity of galactic coronae is low 
primarily because the central density of hot gas has been reduced. 
This, in turn, can be interpreted in terms of a raised adiabat (en- 
tropy) for the gas. It is important to establish when the entropy was 



raised and, if possible, by what mechanism. Fig.|10| shows the me- 
dian temperature (left-hand panel) and entropy (right-hand panel) 
histories of all gas that is in the X-ray luminous phase at z — 0, 
binned by halo mass (also at z — 0). The sharp rise in both quanti- 
ties at z — 9 is caused by the onset of Hi reionisation, which is as- 
sumed to occur instantly and uniformly. In absolute terms, however, 
the entropy and thermal energy injected into the gas from reionisa- 
tion are small compared to the final entropy and temperature of the 
hot gas. However, a second sharp rise in both quantities is evident 
over the epoch 1 < z < 3, corresponding to the peak in A/* for our 
galaxies (which we investigate in greater detail in § |4.3.2fr . Unfortu- 
nately, this rise alone does not enable us to distinguish between star 
formation and feedback as the primary mechanism by which the 
median entropy of hot gas is raised, since feedback is intimately 
tied to star formation, both spatially and temporally (both in real 
galaxies and our simulations). 

In order to obtain insight into which physical mechanism is 
primarily responsible for lowering the central gas density of the 
simulated disc galaxies, we plot in Fig.fTTJthe sorted gas entropy as 
a function of enclosed gas mass. This is calculated for each galaxy 
by ordering the hot gas halo particles by their entropy and plotting 
the entropy, S, against the total mass of gas of particles with en- 
tropy lower than S. In this diagram, cooling selectively removes 
the lowest entropy gas from the hot phase (converting it into cold 
gas and stars), which results in a shift in the S(< M gas ) curves to 
the left (see |Bryan|2000||Voit et al.|2002| >. Feedback, by contrast, 
actually heats the gas (shifts the curves upward) but, if it is ener- 
getic enough, it may also eject gas from the system (shift to the 
left). 

The solid coloured curves in Fig. fTTJ represent the median 
S(< Mg as ) profiles for various mass bins. The solid black curve 
represents the self-similar result obtained from non-radiative cos- 
mological simulations (see |Voit et al.|20 05). To determine whether 
radiative cooling/star formation is the dominant process, we trun- 
cate the self-similar S(< M gas ) distribution (i.e., shift to the left) 
using the calculated stellar fractions from the simulated galaxies 
(dashed curves). This is done by removing all gas with entropies 
lower than the entropy that encloses / gas = /star- (The remaining 
gas mass fraction is then just the original gas mass fraction mi- 
nus the stellar mass fraction.) This procedure yields the expected 
S(< A/g as ) distribution after cooling has removed the lowest en- 
tropy gas from the hot phasaH 

A comparison between the dashed and solid curves demon- 
strates that it is only for the most massive haloes in our sample 
that the truncated self-similar distribution approaches the actual 
S(< A/gas) distribution in our simulations. For these systems, 
therefore, radiative cooling/star formation alone is sufficient to ex- 
plain the high entropy of the halo gas. For lower mass systems, on 
the other hand, selectively removing the lowest entropy gas does 
not by itself raise the entropy of the hot gas sufficiently to re- 
produce the simulated S(< M gas ) distribution. For a halo with 
Af2oo ~ 10 12 Mq, for example, the central entropy (within, say, 
/ g as/(fib/^m) = 0.1) can be raised by a factor of w 5 by radia- 
tive cooling/star formation, but another factor ~ 5 is required to 
match the simulations. Thus, for galaxies with masses similar to 



7 More precisely, this represents an upper limit to the entropy distribution 
after cooling has selectively removed the lowest entropy gas. It represents 
an upper limit as it has been implicitly assumed that any gas that has not 
cooled out of the hot phase flows adiabatically inward to replace the gas 
that has cooled out. In reality, however, some of the remaining hot gas will 
have its entropy reduced by cooling. 
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Figure 10. The median temperature and entropy histories of gas that resides in the hot coronae of our simulated ~ L* disc galaxies at z = 0. Galaxies are 
binned as per Fig. 8 The profiles are normalised to the present-day virial temperature and entropy, respectively, of each system (S200 <* ^200 °^ ^200 )■ ^ ne 
rapid elevation of both quantities at z = 9 signals the onset of Hi reionisation. Both quantities steadily rise to z = 0, experiencing the most rapid increase 
over the interval z ~ 1 — 3, corresponding to the peak in the star formation history of the simulations (see C09). 



the Milky Way and lower, feedback (at z ~ 1 — 3) begins to domi- 
nate entropy production. This is consistent with the reduced overall 
baryon mass fractions of haloes with this mass and lower seen in the 
top-left panel of Fig. [7] Note that feedback from accreting super- 
massive black holes, which is neglected in the present study, would 
presumably reduce the importance of the cooling mechanism for 
the more massive systems. 



4.3 The physical nature of X-ray-emitting gas 

We have demonstrated that star formation (via radiative cooling) 
and, more importantly, feedback from star formation, particularly 
during the epoch corresponding to the peak in the cosmic star for- 
mation rate density (1 < z < 3), determine the state of the high 
entropy, low density X-ray luminous coronae of Milky Way-like 
galaxies in the GIMIC simulations. This contrasts with the picture 
posited by WF91, in which the dense corona that forms around a 
galaxy is assumed to be isothermal and to trace the dark matter dis- 
tribution. In their model a cooling flow is established which fuels 
the ongoing star formation necessary to replenish stellar mass lost 
to the ISM by stellar evolution, and to offset the transfer of mass 
from discs to bulges that results from mergers and bar instabilities. 
As seen in Fig. B] despite having hot gas reservoirs that are less 
dense and have longer cooling times than assumed in the WF91 
model, our sample of galaxies still exhibits ongoing star formation. 
This raises two key questions regarding the coronae of the simu- 
lated galaxies: i) is the hot gas of the corona inflowing (accreting) 
or outflowing (in a wind), and ii) are galaxy discs replenished by 
ongoing star formation fed by the accretion of hot or cold gas? We 
explore these questions in order. 



4.3.1 X-ray luminous gas: inflowing or outflowing? 

Since the median temperature and entropy histories shown in 
Fig. [TO] evolve smoothly, one might naively rule out the hypothe- 
sis that the majority of the gas surrounding our simulated galax- 
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Figure 11. The sorted gas entropy as a function of enclosed gas mass frac- 
tion (solid curves). This is calculated for each galaxy by ordering the hot gas 
halo particles by their entropy and plotting the entropy, S, against the total 
mass of gas of particles with entropy lower than S. The solid black curve 
represents the self-similar result obtained from non-radiative cosmological 
simulations (see Voit et al. 2005 1. The dashed curves represent the expected 
S(< Af gas ) distribution after cooling has removed the lowest entropy gas 
from the hot phase, which is calculated by removing from the self-similar 
distribution all gas with entropy that encloses / gas = / s t a r (see text). The 
similarity of the dashed and solid curves for high halo masses implies that 
radiative cooling/star formation is primarily responsible for 'raising' the en- 
tropy of the (remaining) hot gas. For low mass haloes, however, feedback 
is also required to raise the adiabat. This mode of entropy production dom- 
inates for systems with masses similar to the Milky Way and lower. 
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ies was recently deposited into the hot, X-ray luminous phase by 
supernova-driven winds. To assess this hypothesis more precisely, 
however, we resort to a more direct test of the inflow/outflow sce- 
nario. 

It is possible to measure the fraction of X-ray luminous gas in 
outflows because the simulation code tracks when, if ever, a bary- 
onic particle last attained a density greater than the threshold for 
star formation (hh =0.1 cm -3 ). Particles with density greater 
than this threshold are considered part of the ISM, and are subject 
to an imposed equation of state (EOS), P oc p 7EOS , where our 
choice of 7eos = 4/3 ensures that both the Jeans length and the 
ratio of the SPH kernel to the Jeans length are independent of den- 
sity ( Schaye & Dalla Vecchia 2008 1. In order to transit into the hot 
phase, ISM particles (i.e., those on the EOS) must have formed part 
of an outflow driven by winds. 

The left-hand panel of Fig. [12] shows the mass fraction of 
presently X-ray luminous gas that has been in the ISM, i) at any 
time in the past (black dots), and ii) within the last 2 Gyr (i.e. z < 
0.17, red circles), a period roughly corresponding to the dynami- 
cal time of the haloes in our sample. Typically, only 10 — 20 per 
cent by mass of the X-ray luminous gas presently in galactic coro- 
nae ever reached star forming densities. The fraction that attained 
those densities within a dynamical time of the present day - and 
could potentially therefore form part of an outflow at the present 
epoch - is smaller still, typically less than 2 per cent. However, an 
interesting trend is seen as a function of halo mass, such that lower 
mass haloes do show a large scatter in this quantity. In a small num- 
ber of cases, outflows dominate the overall mass of X-ray luminous 
gas. 

As shown in the right-hand panel, the luminosity fraction con- 
tributed by these particles is disproportionately large for their mass. 
This is not unexpected since much of this gas remains relatively 
dense at the centre of the system even after being heated. More- 
over, having formed part of the ISM, these particles are typically 
more metal rich, boosting their emissivity. Typically, however, gas 
that has recently been heated out of the ISM and into the X-ray 
luminous phase contributes less than half of the overall soft X-ray 
luminosity of galactic coronae. 

A more direct measure to establish the importance of present- 
day outflows to the overall emissivity of hot coronae is to identify 
hot gas particles with positive radial velocities, f ra d, with respect 
to the galactic centre. In Fig.[T3]we show the fractional contribution 
of these particles to i) the X-ray luminous gas mass of haloes (left- 
hand panel) and ii) the soft X-ray luminosity of haloes. We show 
the contribution from particles in winds that exceed three velocity 
cuts, measured in terms of the halo circular velocity at the virial ra- 
dius: «2oo/10 (black dots), «20o/3 (red dots) and «2oo (blue dots). 
Only a modest fraction of the hot gas mass exhibits an apprecia- 
ble positive radial velocity, and very little approaches the circular 
velocity. It is therefore safe to conclude that all hot gas in the sim- 
ulated systems remains gravitationally bound to the halo, since the 
escape velocity is typically a factor of a few greater than the circu- 
lar velocity. As also concluded from inspection of Fig. [12] most of 
the mass of X-ray luminous gas is inflowing (or near hydrostatic), 
albeit with a trend in halo mass such that outflows contribute more 
to the X-ray luminous mass fraction in less massive haloes. Again, 
in common with Fig. [12] outflows contribute a disproportionate 
fraction of the overall X-ray luminosity because their density and 
metallicity are typically greater than for inflowing gas. 

In §[5] we demonstrated that the GIMIC simulations reproduce 
the observed correlation between the present-day star formation 
rate, M*, and the soft X-ray luminosity, I/x,o.5-2.okcV (see Fig.pl. 



This observed scaling has been interpreted by several authors as 
evidence that the formation of hot coronae is linked to supernova- 
driven feedback associated with star formation. However, we have 
just seen that typically less than half of the X-ray luminosity is 
contributed by recently heated, outflowing gas. Moreover, no tran- 
sition in this trend is seen close to «2oo ~ 250 km s" 1 (denoted 
by a down arrow in Fig. |13[ l, the scale at which haloes exhibit a 
sharp transition in star formation efficiency (see Fig. 8 and associ- 
ated discussion in C09). The absence of such a trend indicates that 
present-day star formation (and its associated supernova feedback) 
in general does not drive the X-ray luminosity. 

The M t — Lx,o.5-2.okcV scaling has a simple, alternative ex- 
planation: both quantities scale with the virial mass of the halo, 
M2oo- A direct illustration is provided by Fig. [14] in which X-ray 
luminosity is plotted as a function of the halo specific star forma- 
tion rate of the galaxy, M±/M2oo- Note that, in this case, we plot 
the bolometric X-ray luminosity 8 ] as computed by APEC, to pre- 
vent a spurious drop in Lx for cases where a significant fraction of 
the emission falls outside of the soft X-ray waveband. With the halo 
mass dependence factored out of the star formation rate, it becomes 
clear that in the most luminous systems (Lx.tot > 10 41 erg s _1 ), 
the X-ray luminosity is weakly anti-correlated with M*/M2oo- For 
less luminous systems a modest positive correlation is visible, al- 
beit with considerable scatter. 

These trends are consistent with the conclusion we drew from 
Figs. [12] and [T3] namely that in massive systems the X-ray lumi- 
nosity is dominated by quasi-hydrostatic or inflowing gas, whilst 
in smaller counterparts outflowing gas contributes a significant, but 
not dominant, fraction of the luminosity. It is also noteworthy that 
observations are necessarily luminosity weighted. As a result, X- 
ray data yield a rather unrepresentative view of the hot gas dis- 
tribution, since the (? dependence of Lx leads to a strong bias- 
ing of the luminosity towards the halo centre, relative to the over- 
all mass distribution of X-ray luminous gas. We find that, inde- 
pendently of halo mass, half of the X-ray luminosity is contained 
within ~ 0.1 — 0.2r2oo, whilst half of the X-ray luminous mass 
is contained within ~ 0.6 — 0.7r2oo- For this reason, the emis- 
sion detected by X-ray telescopes is disproportionately weighted 
towards outflowing gas. Even with this, we predict that at least half 
of the X-ray luminosity detected around disc galaxies is from quasi- 
hydrostatic inflowing gas. High resolution spectroscopic X-ray ob- 
servations enabled by future telescopes, such as NeXT/ASTRO-H 
(Takahashi et al. 2008) and eventually the International X-ray Ob- 
servatory (IXO), have the potential to map directly the velocity 
structure of this gas, thus providing constraints on the fraction of 
hot gas that is entrained in outflows. 



4.3.2 What fuels star formation in discs: hot vs cold accretion? 

In Fig. [8] we showed that, at z = 0, the cooling time of hot gas 
surrounding L* disc galaxies is considerably longer than predicted 
by WF91. This deviation from the analytic model occurs because, 
in the hydrodynamic simualtions, star formation and supernova- 
driven feedback at z ~ 1 — 3 raise the adiabat of hot gas prevent- 
ing its compression to high densities at later times. Since we find 
significant differences between the simulations and the canonical 
analytic picture of WF91, we are motivated to investigate whether 



8 Whilst semantically inaccurate, the term 'bolometric X-ray luminosity' is 
frequently used in the field. We adopt it here to mean the luminosity yielded 
by a plasma emission model, integrated over the 0.02-50 keV energy range. 
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Figure 12. The fraction of X-ray luminous mass (left-hand panel) and soft X-ray luminosity (right-hand panel) contributed by gas in our simulated galaxies 
that has previously been in the star-forming 'equation of state' (EOS) phase that mimics the thermodynamics of the ISM. By definition, transit from the ISM 
to the hot phase is via an outflow. Each galaxy contributes one black dot and one red circle to each plot; black dots denote the fractions due to gas that has 
been in the ISM at any time in the past, whilst red circles show the fractions for gas that has been in the ISM in the last 2 Gyr, corresponding to a dynamical 
time for L* galaxies. Above M200 ~ 10 12 Mq, only 10 — 20% of the hot gas mass has ever been processed through the ISM, and only a few per cent has 
been within the last 2 Gyr. However, this gas contributes a disproportionately high fraction of the coronal X-ray luminosity, owing to its centrally concentrated 
distribution and relatively high metallicity. 
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Figure 13. The fraction of X-ray luminous mass (left-hand panel) and soft X-ray luminosity (right-hand panel) contributed by gas in our simulated galaxies 
that has a positive radial velocity at z = 0. Each galaxy contributes three dots (one each of black, red and blue) to each plot. We delimit 'slow' (black dots), 
'intermediate' (red dots) and 'fast' (blue dots) outflows with velocity cuts scaled to the system circular velocity, as shown in the legend. The down arrow at 
"200 = 250 km s _1 marks the approximate scale above which supernova-driven winds become inefficient at regulating star formation (see C09). The plots 
are consistent with Fig. |12| relatively little X-ray luminous mass is found in fast outflows. However, outflowing gas is typically more luminous (per unit mass) 
than hydrostatic or inflowing gas, owing to its centrally concentrated distribution and relatively high metallicity. The absence of a transition in behavour at 
v c — 250 km s~ 1 indicates that the X-ray luminosity of the simulated galaxies is not in general driven by the present-day star formation rate. 



the cooling of gas from a hot, quasi-hydrostatic reservoir is really 
what fuels star formation in galaxy discs at late times, as posited 
by WF91. As shown in Fig.B] our simulated galaxies exhibit star 
formation rates that are consistent with observations; this indicates 
that the gas consumed by star formation in the simulations is ef- 



ficiently replenished, or otherwise the cold gas in discs would be 
rapidly consumed. 

It is less clear how this replenishment actually takes place. 
Recent studies employing hydrodynamic simulations have claimed 
that a significant fraction of a galaxy's stellar mass forms from gas 
that is channelled to the halo centre without ever experiencing a 
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Figure 14. The 'bolometric' X-ray luminosity (0.02-50 keV) of hot gas 
coronae surrounding our simulated galaxies, as a function of their specific 
star formation rate (sSFR; M*/Af2no)- Although we saw in Fig.plthat X- 
ray luminosity tracks the star formation rate, both quantities correlate with 
the system virial mass, M200 ■ With this scaling factored out, a weak corre- 
lation exhibiting a large scatter is evident for less massive systems, whilst 
for more massive systems the X-ray luminosity weakly anti-correlates with 
the sSFR. The Lx — M* correlation is therefore driven primarily by the 
scaling of both quantities with M200 ■ 



sustained accretion shock (Keres et al. 2005] |Dekel et ar|2009| >. 
The thermalisation of the potential energy of this material occurs 
at significantly greater densities than those near the virial radius, 
drastically shortening its cooling time and suppressing the forma- 
tion of a hot coronal phase. This scenario has spawned the term 
'cold accretion'. It corresponds roughly to the regime identified by 
WF91 in which the cooling time of the gas is much shorter than the 
dynamical time of the halo (see Fig. 2 of WF91) although WF91 
considered only a spherically symmetric case whereas much of the 
'cold accretion' in the simulations seems to be associated with fil- 
aments. 

It is instructive to assess, for our sample of simulated galax- 
ies, the fraction of star formation from gas that has arrived at the 
disc i) having shocked and subsequently cooled (as posited by 
WF91 and similar models - the 'hot mode') and ii) without hav- 
ing experienced a sustained accretion shock ('cold mode'). This is 
relatively straightforward to determine since our simulation code 
tags star particles with both their formation redshift and the max- 
imum temperature attained by the gas particle from which they 
formed. We therefore adopt a simple definition of the two modes: 
stars that formed from gas whose temperature never rose above 
T = 2.0 x 10 5 K (the minimum temperature we consider for X- 
ray emission) are considered to have formed via the cold mode; all 
others are considered to have formed via the hot mode. Note that 
this simple definition of 'cold' and 'hot' is not identical to that of 
e.g. |Keres et al.| ((2005) and tends to exaggerate the amount of cold 
accretion since at least some of the gas below the temperature cut 
is likely to have been shock-heated to the virial temperature. 

The left-hand panel of Fig.[l5]compares the instantaneous star 
formation rate fuelled by the hot (solid lines) and cold (dotted lines) 
modes, as a function of redshift for all progenitors of our simulated 
galaxies. The galaxies are binned according to their virial mass at 



z = and their median star formation histories are plotted. The 
cold mode dominates early star formation in all mass bins, but in 
all cases the hot mode overtakes the cold mode during the interval 
1 < z < 2 (which corresponds roughly to half the age of the Uni- 
verse). The right-hand panel of Fig.[l5]shows the stellar mass frac- 
tion formed from gas accreted in each mode up to a given redshift. 
The fraction of stars formed in each mode is almost independent of 
halo mass, in the range considered here. In all cases, most of the 
stars present in these haloes today formed from gas accreted in the 
hot mode. The fraction ranges from ~ 65% for the least massive 
haloes plotted in the figure to ~ 55% for the most massive. Weight- 
ing by the relative numbers of haloes of different mass, the fraction 
of stars formed from hot mode accretion in haloes with mass above 
the minimum shown in the figure (M200 = 10 Mq) is ~ 54%. 
The values of these quantities no doubt depend, to some extent, on 
the particular feedback scheme adopted in the simulations, but any 
scheme that gives a reasonable galaxy stellar mass function is likely 
to produce similar numbers. The dependence of various modes of 
star formation on the feedback model is being investigated in the 
OWLS simulations by van de Voort et al. (in preparation). Beyond 
this, it will be important to check that the results are independent 
of the numerical hydrodynamics technique employed, whether it be 
SPH, AMR or a moving mesh scheme. 

The dominance of hot mode star formation at low redshift is 
of particular relevance here because it supports the simple assump- 
tion of the WF91 model that disc formation at the present epoch 
is fuelled by the cooling and condensation of gas from a hot halo. 
This is in spite of the fact that, as we have shown here, the gas in 
the hot haloes lies on a higher adiabat than assumed by WF91 (i.e. 
it is hotter, less dense, and has a longer cooling time). Note that this 
general picture is not inconsistent with the observation that extra- 
planar high-velocity clouds (HVCs) are abundant in the Milky Way 
and nearby disc galaxies and may be transporting a significant mass 
of cold gas to their discs (e.g.|Oort|1970|[Kerp et al.|1996||Wakker 



& van Woerden 1997 ; Wakker et al. 1999 ; Tripp et al. 2003 . 



Thilker 



et al. 2004; Miller et al. 2009). Such clouds may well have formed 



from thermal instability and fragmentation of hot coronal gas (e.g. 
|Maller & Bullock|2004| [Booth & Theuns|2007| >. 

4.4 Comparison with other theoretical work 

At this juncture it is instructive to draw comparisons with previ- 
ous attempts to compute the X-ray luminosity of hot coronae of 
disc galaxies at z = 0. The simulations discussed in this context 
by Rasmussen et al. (2009), initially presented by Sommer-Larsen 
|et al.| i 2005 1 and Romeo et al.| (2006), and which were an updated 
version of those discussed in a similar context by |Toft et al.| (2002 ), 
arrived at seemingly similar conclusions to those presented here. 

However, the similarity between the two sets of conclusions 
is probably coincidental because of the very different treatments 
of baryonic physics in the two studies. In particular, the super- 
nova feedback scheme implemented by Rasmussen et al. (2009) 
was only effective at driving galactic scale outflows at z > 4 — 5. 
As a result, star formation in their galaxies was not self-regulating 
and ended up with stellar mass fractions (/* = M+/M200) roughly 
twice as large as those in our GIMIC sample (J. Sommer-Larsen, 
private communication). Such high stellar fractions are ruled out 
by analyses of the Sloan survey (e.g. |Guo et al.| 2010), indicating 
that the X-ray luminosities produced by galaxies in the simulations 
of |Toft et ak1l2002| ) and |Rasmussen et al.| ( [2009^ are lower than 
predicted by WF91 primarily because of overcooling and exces- 
sive star formation. According to our results, these studies under- 



Age of the Universe [Gyr] 



X-ray coronae in simulations of disc galaxy formation 17 



Age of the Universe [Gyr] 




2.1 1.5 1.2 0.9 0.8 
I I I I r 



"g 0.8 



[12.50,12.75] 
[12.25,12.50] 
[12.00,12.25] 
[11.75,12.00] 
[11.50,1 1.75] 



Hot 
Cold 



,... s ;;;::: : "' 




Figure 15. The star formation history of all progenitors of our simulated galaxies, binned as in Figs.[8]and |10| decomposed into hot (solid lines) and cold 
(dotted lines) modes (see text for details). The left-hand panel shows the instantaneous star formation rate in each mode, whilst the right-hand panel shows 
the fraction of the total stellar mass in each bin that has formed in each mode by the given redshift. In every mass bin, star formation fuelled by cold mode 
accretion dominates at early epochs, whilst star formation from the hot mode becomes dominant at later epochs. The increase in hot mode star formation rate 
occurs sufficiently early that in every mass bin the stellar mass formed by the hot mode overtakes the cold mode mass at z ~ 1. Most stars today in our sample 
of ~ L+ galaxies formed from gas accreted in the hot mode. 



estimate the role of feedback, which we found to be dominant, in 
explaining why the X-ray luminosities of disc galaxies are much 
smaller than predicted by WF91. 

In principle, the effects of radiative cooling and feedback can 
be included in semi-analytic models, as shown by Bower et al. 
(2008), who attempted to construct a self-consistent model of the 
heating of the intracluster medium by AGN and its effect on the 
subsequent cooling of the gas. This study focussed on the X-ray 
properties of groups and clusters, but the methodology is applica- 
ble to galaxies as well. 



5 SUMMARY AND DISCUSSION 

The existence of hot, X-ray luminous gas reservoirs surround- 
ing present day Milky Way-like galaxies is a key prediction of 
the galaxy formation model developed by White & Frenk (1991 
WF91). This prediction is important because the WF91 model, and 
the many analytic and semi-analytic models based on elements 
of it, successfully reproduce a broad range of observed galaxy 
properties. However, evidence in support of this particular predic- 
tion has, in general, not been forthcoming. In this paper, we have 
used the Galaxies-Intergalactic Medium Interaction Calculation 
( |Crain et al.||2009| C09) to investigate this challenge to the stan- 
dard picture of galaxy formation using altogether different tech- 
niques. About half the galaxies that form in these five simulations 
are disc-dominated (an interesting fact in its own right which we 
will explore in detail in a further paper). Of these, 458 have stellar 
masses greater than 10 10 hT 1 Mq and make up the sample we have 
investigated here. 

Benson et al. (2000 BOO) searched for the predicted extended 
emission around three promising candidate disc galaxies using the 
ROSAT telescope, but they failed to detect any. Subsequent exten- 
sive surveys with the XMM-Newton and Chandra telescopes finally 
succeeded in detecting diffuse emission in the soft (0.5-2.0 keV) X- 



ray band around some disc galaxies, but at much lower levels than 
predicted by WF9 1 . The larger number of non-detections in these 
surveys has tightly constrained upper limits on this kind of emis- 
sion. The observational data now offer a clearer picture than was 
available to BOO who were the first to highlight the potential impli- 
cations of the absence of X-ray emission for galaxy formation the- 
ory. In sum, the data indicate that the soft X-ray luminosities of disc 
galaxies are one to two orders of magnitude lower than predicted 
and that they broadly correlate with both the i\~-band luminosity 
and the star formation rate of the galaxy. 

The GIMIC hydrodynamic simulations make different assump- 
tions and approximations to those of the WF91 model. However, 
we find that these simulations also generate X-ray coronae around 
disc galaxies, although these are much weaker that predicted by 
WF9 1 for the reasons summarized below. Many of the bulk X-ray 
properties of the simulated galaxies agree with observations. In par- 
ticular, the 458 galaxies in our sample exhibit similar scaling rela- 
tions to the data. For this comparison, we used a compilation of 



X-ray measurements and upper limits from several studies iBen 



son et al.|200"0)|S trickland e t al.|2004|fW ang 2005; Tullma nn et al. 
2006||Sun et al.|2007||Owen & Warwick|2009| |Rasmussen et ~ 



2009 ) and extracted i\"-band luminosites from the 2MASS database 
and disc rotation velocities (derived from inclination-corrected 21- 
cm measurements) from the HyperLeda database. We find that 
the galaxies in the simulations broadly reproduce the observed scal- 
ing and scatter of the Lx — Lk, Lx — «rot and Lx — M± relations 
at 2 = 0. 

There are two main reasons why the X-ray emission from the 
GIMIC galaxies is much weaker that predicted by WF91. Firstly, 
the mass of hot gas within the galaxies' haloes is, typically, lower 
than in the WF91 model. This is partly because some of the gas has 
been consumed into stars and partly because some gas has been 
blown out in winds. Secondly, and most importantly, the distribu- 
tion of the hot gas is much less concentrated than the distribution 
of dark matter, violating one of the main assumptions of the WF91 
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model. The gas is more extended than expected because its entropy 
has been significantly raised. For galaxies like the Milky Way and 
fainter, the increase in entropy is mostly due to the injection of en- 
ergy generated by supernovae. The gas is placed in a high adiabat 
at early times, z ~ 1 — 3, when the star formation peaks. 

The entropy of the corona can also be affected by radiative 
cooling. Cooling acts to reduce the gas temperature by selectively 
removing the lowest entropy gas and locking it into stars, thus rais- 
ing the mean entropy of the remaining coronal gas (e.g. Voit et al. 
2002). The effect of radiative cooling in this respect is therefore 
not dissimilar to that of outflows. However, we have shown that the 
reduction of the hot gas density by radiative cooling is only signif- 
icant for the most massive galaxies in our sample. In this regime, 
however, our simulations overestimate the effects of radiative cool- 
ing because they do not include additional energetic feedback pro- 
cesses such as AGN. This shortcoming is manifest, for example, in 
the high stellar fractions of the largest galaxies which are about a 
factor of 2 larger than inferred observationally (Guoetal.|2010| l. By 
contrast, the median stellar fractions of the Milky Way-like galax- 
ies and smaller in our simulations are in very good agreement with 
these data, although the scatter in the simulations appears perhaps 
to be too large. Simulations that include the effects of AGN are re- 
quired to understand the X-ray properties of the gas associated with 
the most massive galaxies. 

The main difference between the hydrodynamic simulations 
and the semi-analytic models based on the WF91 formalism is that 
the simulations follow the kinematics and dynamics of the gas self- 
consistently. Although outflows are often included in semi-analytic 
models (see e.g.[K auffman n~eTaT.|1999bHCole et al.|20 00; Bower 
et~aLl|200T| [Baugh et al.||2005| |Bower et al.||2006||Croton et al.| 



2006b. Bertone et al. 2007), the readjustment of the gas to energy 



injection is not generally taken into account. One exception is the 
recent work by [Bower et al. (2008) which attempts to follow self- 
consistently the heating of the intracluster medium by AGN and its 
effects on the subsequent cooling of the gas. The key contribution 
of the full hydrodynamic treatment is the ability to follow the way 
in which the entire structure of the gas corona adjusts in order to 
accommodate the increase in entropy associated with the energy 
injection by supernovae (or AGN). It is important to note, how- 
ever, that the GIMIC simulations assume a specific model for energy 
feedback. Although this appears well motivated, it is essential to 
investigate how the X-ray properties of the gas depend on the im- 
plementation of energy feedback processes. Whilst methods such 
as that developed by Bow er et al.| 12008 ) provide a blueprint for 
starting to explore the reaction of gas to energy injection in a com- 
putationally inexpensive fashion, we note that model-independent 
conclusions can only be arrived at with suites of hydrodynamical 
simulations featuring a wide range of possible feedback treatments 
(an approach adopted by the OWLS simulations, albeit at lower res- 
olution than required here). We intend to pursue this methodology 
with high-resolution resimulations of individual galaxies in future 
studies. 

In spite of the approximations inherent in the analytic treat- 
ment of gas dynamics in the WF91 model, key results of this 
model are borne out by the GIMIC simulations. The predicted quasi- 
hydrostatic hot coronae around L* disc galaxies do indeed form, as 
anticipated, by the shock heating of infalling gas as the dark matter 
halo develops. In all but a few low mass galaxies, emission from 
these coronae dominates the galactic diffuse X-ray luminosity. 

Galactic winds contribute to the X-ray luminosity but, in gen- 
eral, their emission is subdominant. Winds represent a small frac- 
tion of the coronal gas, but since they tend to be hotter, denser and 



more metal-rich than the quasi-hydrostatic component, they con- 
tribute disproportionally to the X-ray luminosity. Since most X- 
ray observations target the central parts of disc galaxies, where the 
wind gas is concentrated, and since the observed X-ray luminos- 
ity correlates well with the star formation rate, it is easy to misin- 
terpret the importance of winds to the overall X-ray emission. In 
fact, in our simulations the X-ray luminosity also correlates with 
the star formation rate, just as in the data, even though the emis- 
sion is dominated by the quasi-hydrostatic corona. This correlation 
occurs because both the star formation rate and the total amount of 
hot gas depend on the mass of the halo. The simulations do pro- 
duce a few examples of bright galaxies in which outflowing gas 
contributes the bulk of the X-ray luminosity. These may be coun- 
terparts of starburst galaxies such as M82 (whose intense star for- 
mation seems to be associated with interaction with M81, sec e.g. 
Chynoweth et al.|2008 ). In our simulations, star formation in L* 
disc galaxies today is mostly fuelled, as found by WF91, by the 
cooling of gas out of the hot corona. This can occur either through 
a smooth cooling flow (e.g. |Thomas et al.|1986) or through frag- 
mentation into high-velocity clouds (e.g. Mailer & Bulloc kf2004[ 
|Booth & Theuns|2007) . 

Our simulations indicate that hot outflows do not, in gen- 
eral, dominate the overall soft X-ray luminosity of isolated disc 
galaxies. We therefore expect gravitational inflow to dominate over 
supernova-driven outflows, and the majority of a galaxy's soft X- 
ray luminosity to stem from hot, quasi-hydrostatic gas at low sur- 
face brightnesses. This physical picture is, in principle, testable 
with future X-ray observatories with instrumentation of high sen- 
sitivity and spectral resolution, such as NeXT/ASTRO-H and 1XO, 
because coronal gas dynamics follow directly from X-ray line di- 
agnostics. We therefore expect these facilities directly to verify, or 
rule out, the key role in the formation of disc galaxies that theoret- 
ical models assign to hot gaseous coronae. 
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Al. The correspondence between the soft X-ray luminosity 
-2.0 kcv) °f galaxies in the intermediate-resolution — 2cr GIMIC 



Figure 

(-^X.0.5 

simulation, and that of their counterparts in the high-resolution realisation 
of the same region. The dashed line denotes the 1 : 1 relation, and the dotted 
lines are offset higher and lower by a factor of 2. 
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APPENDIX A: CONVERGENCE 

Our sample of simulated galaxies is drawn from the five 
intermediate-resolution GIMIC simulations, which form a complete 
set at z — 0. In addition, the GIMIC suite includes one simula- 
tion to z — 0, that of the low density — 2a region, with 8 times 
better resolution. This simulation allows an assessment of the nu- 
merical convergence of our results. In C09 we used this and other 
high-resolution simulations at high redshift to show that the star 
formation rate density in the intermediate-resolution simulations 
has converged for z < 6, and that the halo specific star forma- 
tion rate at z — 0, M+/M200, has converged for haloes with 
W200 > 100 km s" 1 . 

The main quantity of interest here is the soft X-ray luminosity 
of galaxies at z — 0. For each galaxy in our sample drawn from the 
—2cr region at intermediate resolution, we have identified its coun- 
terpart in the high-resolution realisation. We assess convergence by 
comparing the soft X-ray luminosity of the objects at intermediate 
and high resolution in Fig. |Al| The dashed line traces the locus of 
L™ t_reB = L x 1_rcs , whilst the upper and lower dotted lines are 
offset higher and lower by a factor of 2. 

The data exhibit excellent convergence properties, with more 
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than 50 percent of systems having soft X-ray luminosities in the 
intermediate- and high-resolution runs that agree to within a factor 
of 2, and 85 percent that agree to within a factor of 5, over 4 decades 
in Lx- The degree of convergence in Lx is more than adequate 
for the purposes of this work. For example, the spread in Lx for 
galaxies of a given Lk in the observational data is typically much 
greater than the typical difference in Lx for a given galaxy sim- 
ulated at intermediate- and high-resolution. Thus, we expect that 
our main conclusions would remain unchanged if they had been 
obtained from a full set of 5 GIMIC regions at high resolution. 



